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The transition-matrix (T-matrix) approach provides a general formalism to study scattering prob- 
lems in various areas of physics, including acoustics (scalar fields) and electromagnetics (vector 
fields), and is related to the theory of the scattering matrix (S'-matrix) used in quantum mechanics 
and quantum field theory. Focusing on electromagnetic scattering, we highlight an alternative for- 
mulation of the T-matrix approach, based on the use of the reactance matrix or A'-matrix, which 
is more suited to formal studies of energy conservation constraints (such as the optical theorem). 
We show in particular that electrostatics or quasi-static approximations can be corrected within 
this framework to satisfy the energy conservation constraints associated with radiation. A general 
formula for such a radiative correction is explicitly obtained, and empirical expressions proposed in 
earlier studies are shown to be special cases of this general formula. This work therefore provides a 
justification of the empirical radiative correction to the dipolar polarizability and a generalization 
of this correction to any types of point or body scatterers of arbitrary shapes, including higher 
multipolar orders. 

PACS numbers: 42.25.Fx, 78.20.Bh, 41.20.-q, 11.55.-m 



I. INTRODUCTION 

Radiative reaction, also known as radiation damping, 
refers to the fact that the electromagnetic (EM) field cre- 
ated by a charge or emitter must affect its own dynamics 
(e.g. motion or power radiated) pp. When applied to ele- 
mentary charges [l2j , no satisfactory classical treatment of 
this effect has been found [3J, yet the radiative reaction 
is at the core of the concepts of self-energy and renor- 
malization in quantum electrodynamics (QED) [1]. In 
fact, using a Green's function approach and regulariza- 
tion techniques akin to those of QED, a classical treat- 
ment of the radiative reaction for point electric dipole 
scatterers can be obtained [4]. 

Interestingly, an equivalent result had been obtained 
heuristically by adding a reaction field postulated from 
simple energy conservation arguments [5]. These argu- 
ments were inspired by research into simple models of the 
optical properties of sub-wavelength particles, notably 
for applications in plasmonics and surface-enhanced Ra- 
man spectroscopy [5J, and we briefly present a similar 
derivation here. The main idea is to use the solution 
of the electrostatics problem to derive an approximate 
dipolar polarizability ao for the particle, ao, assumed 
isotropic here for simplicity, defines the electrostatics re- 
sponse of the particle and is such that a uniform ex- 
ternal electrostatic field Eo induces a dipole moment 
Po = aoEo- In such an electrostatics problem, the 
power absorbed by the particle equals the work done 
by the external field on the charges [3J and is therefore 
P a ° bs - (l/2Vlm(a )|E | 2 . 

In the electrostatics, or quasi-static, approximation, 
also often called Rayleigh approximation (see e.g. Chap. 
5 in Rcf. [7]), the far- field optical response of a sub- 
wavelength scatterer to an incident electric field E; nc os- 



cillating at frequency uj is approximated as that of an 
oscillating induced dipole given by p = aEj nc . Note that 
we use complex notations with the exp(— iut) conven- 
tion and also that SI units are used throughout. We 
also define the wave- vector in the medium k\ — y/e~iuj/c. 
The electrostatics approximation consists in approximat- 
ing a by ao, even if the new polarizability a is in prin- 
ciple different from ao because the electrostatics solu- 
tion does not account for radiation. The power ra- 
diated (or scattered) in this approximation is there- 
fore P sca = (wfcJ|a | 2 |E inc | 2 )/(127re ei). The power ab- 
sorbed by the particle is approximated by its electro- 
statics value P a bs = Pabs- The extinguished power is 
the power extracted by this point dipole from the in- 
cident EM field (the work of the field on the dipole) 
and is P cx t = (l/2)odm(a)|Ej nc | 2 in the general case. 
It therefore reduces in the electrostatics approximation 
(a ~ ao) to P ex t = Pabs, which appears to contradict the 
energy conservation condition P cxt = P sca + P a b s - This 
is not so surprising since the electrostatics solution does 
not account for radiation (scattering) effects. In fact, 
there is no contradiction in the strict range of validity 
of the electrostatics approximation, i.e. in the limit of 
vanishing size, as we then have P sca <C Pext, Pabs since 
ao scales with particle volume. Nevertheless, it is use- 
ful in many instances to correct this problem to extend 
the range of applicability of the electrostatics approxi- 
mation. This can be achieved, as proposed in Ref. \S\ 
by defining a radiation-corrected polarizability, which by 
construction enforces the energy conservation condition 

Pext — Psca ~t~ Pabsi i.e.. 

_ . 1 . Im(a RC ) _ . 1 . k\ 

Im rc = 1 RC12 =Im ) + ~^— 1 

a KO |a KU p a O7re ei 

This condition on the imaginary part only is not suf- 
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ficient to define a RC uniquely (unless Kramers-Kronig 
relations [8] are used) and the additional condition that 
Re(l/a RC ) = Re(l/a ) is usually assumed without fur- 
ther justification to obtain the radiative correction to the 
polarizability as: 



V RC 



1 kf 

i . 

a 67re ei 



or equivalently [5J 



„RC 
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(2) 
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This expression can also be derived rigorously in the spe- 
cial case of spherical particles by expansion of the Mie co- 
efficients [7J OH and can be generalized to spheroidal par- 
ticles 03 [10] . This corrected polarizability has been used 
in numerous contexts, including for example surface- 
enhanced Raman scattering 03 [3] , plasmonics [TDHT^] , or 
the discrete dipole approximation [T5HTB] . More recently, 
this radiative correction has also been generalized, using 
again a heuristic approach based on the optical theorem, 
to the case of point magnetic dipole and dipolar scatter- 
ers with magneto-electric coupling [17] . for the study of 
metamaterials. 

In this paper, we propose a general framework for the 
study and further understanding of the concept of ra- 
diative correction, based on a simple reformulation of 
the T-matrix approach to EM scattering [TH] H9] . The 
T-matrix formalism, often used in conjunction with the 
Extended Boundary Condition Method (EBCM), or null- 
field method, was introduced more than 40 years ago 
[20] and is arguably one of the most elegant and efficient 
methods to solve problems of electromagnetic scattering 
by particles of arbitrary shape and size fit?! I2TH27S] . It 
has been applied, for example, to the study of scattering 
by aerosols [27], metallic nanoparticles [28H3T]. and col- 
lections of spheres [32] , and also to more formal studies 
of EM scattering [33]. It has also been used extensively 
for acoustic scattering (26] [34] . 

Our reformulation emphasizes the important role of 
the reactance matrix, or iVmatrix 18, 35] in relation to 
energy conservation and radiative correction. Although 
the if-matrix has been used occasionally in the past in 
the context of the quantum theory of scattering [T5] 1351 - 
[37] , it seldom appears in EM theory. We show that all the 
aforementioned results for the radiative correction in EM 
scattering are special cases of a general formula derived 
in this work. In addition to highlighting the importance 
of the if-matrix for general scattering theory, this work 
therefore provides a formal justification of existing radia- 
tive correction formulae and a generalization applicable 
to any type of point scatterer or particle of arbitrary 
shape. The latter point is a direct consequence of the 
fact that the T-matrix formulation of EM scattering is ex- 
tremely general. It applies to particles of arbitrary shape 
and may also cover for example JTHJ [25J optically-active 
or anisotropic materials, layered particles, and multiple 



scattering by collections of particles. The proposed K- 
matrix reformulation and associated radiative correction 
are therefore applicable to all the aforementioned cases. 

The paper is organized as follows: in Sec. [IT] we briefly 
review the general principles of the T-matrix approach 
to EM scattering. We then introduce in Sec. |III| an alter- 
native, but closely related, formulation of the problem in 
terms of the if-matrix. In Sec. |IV[ we discuss the impli- 
cations of the if -matrix formulation with regard to ra- 
diative corrections and obtain a general formula (Eq. 21 ) 



for the radiative correction in EM scattering. Finally, in 
Sec. [V] we show explicitly how this formula applies to 
specific cases of radiative correction that have been pre- 
sented in the literature, therefore justifying, and in some 
cases extending, these previously empirical results. 



II. THE T-MATRIX APPROACH 
A. Definition of the T-matrix 

We consider the general problem of electromag- 
netic scattering by a body characterized by a lin- 
ear local isotropic relative dielectric function e 2 (possi- 
bly frequency-dependent) embedded in a non-absorbing 
medium of refractive index n\ (and relative dielectric 
function e\ = n\). 

Within the T-matrix approach [HI [22], the EM field 
solution is expanded in a basis of vector spherical wave 
functions (VSWFs) in a similar fashion as for Mie the- 
ory [7]. We here follow the conventions of Mishchenko 
[T5] for the definition of the VSWFs (see Appendix A] for 
details). The incident field Ej, lc and internal field Ej nt 
(the field in the region inside the particle) are regular at 
r = and can therefore be expressed in terms of regu- 
lar VSWFs denoted M^, nJ> 1} . The scattered field E sca 
must satisfy the Sommcrfcld radiation condition and is 
therefore expanded in terms of outgoing spherical waves 

Explicitly: 



VSWFs denoted M^ 3) ,Ni 3) 



Ei„c(r) = E,^Mi 1) (fc 1 r) + 6 y NW(A ;i r), 
Esca(r) = £*lfcMk 3) (fcir) + q u N^(k 1 r), 
Eint(r) = E,c,Ml 1) (fc 2 r) + dM 1} (k 2 r), (4) 



where ki = (2tt /\)^fel (i = 1,2) are the wavevector am- 
plitudes in region 1 (outside) and 2 (inside) and A is the 
excitation wavelength. These expansions can be repre- 
sented as vectors, for example {pu,q v ) = (p, q) for the 
scattered field, where the index v = (n, m) combines the 
total (n) and projected (|to| < n) angular momentum 
indices. The expansion of the incident field (a^, b v ) for a 
given scattering problem is known, with explicit expres- 
sions existing for example for plane waves [19] . 

By linearity of Maxwell's equations, the coefficients of 
the scattered field are linearly related to those of the 
incident field. This can be expressed explicitly by intro- 
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ducing the T-matrix: 



Pi_ Tf . 



(5) 




where T is an infinite square matrix, which can be writ- 
ten in block notation as: 



(6) 



In principle, from a knowledge of the T-matrix (at a 
given wavelength), one can infer the scattering proper- 
ties for any incident excitation. The T-matrix approach 
is therefore particularly suited for computations of the 
scattering properties of a collection of randomly oriented 
scatterers |38) . which is indeed one of the important ap- 
plications of this formalism 19J. 

We note that linear relationships involving the expan- 
sion coefficients of the internal field can also be written: 



and 



Q 



The T-matrix can therefore also be obtained from 



PQ 



(7) 



(8) 



This expression provides the basis for one of the most 
common approaches to calculating the T-matrix in prac- 
tice, namely the Extended Boundary Condition Method 
(EBCM) or Null-Field Method [H [2TJ HH HE] . Within 
this approach, the matrix elements of P and Q are ob- 
tained analytically as surface integrals over the particle 
surface of VSWF cross-products. 



B. Symmetry, unitarity, and energy conservation 

The T-matrix satisfies [23 [H] symmetry relations re- 
lated to optical reciprocity, along with unitarity relations 
related energy conservation, i.e. the fact that the extinc- 
tion cross-section <r ext is the sum of scattering cr sca and 
absorption er a b s (note that this is related to the optical 
theorem [33 HD]). The optical reciprocity relations 
are typically easy to check and enforce as they are re- 
lated (see App. [d]) to ensuring the symmetry of certain 
matrices |21j . The energy conservation condition is in 
general more problematic. It is typically expressed by 
introducing the S'-matrix (scattering matrix) defined as 
S = I + 2T. For lossless (non-absorbing) scatterers (for 
which Im(e2) = 0), it can then be shown that energy 
conservation is equivalent to S being unitary [19, 21 . In 
terms of the T-matrix itself, this results in the somewhat 
more cumbersome condition: 



T + T t _ _2T f T, 



(9) 



which can be viewed as the matrix form of the generalized 
optical theorem fTS] [35] . 



In EM scattering, absorbing or conducting scatterers, 
for which Im(e 2 ) > (we exclude the special case of per- 
fect conductors here), are also often considered and the 
equality above no longer holds. In this general case, the 
inequality <r cxt > <7 sca then requires thatl-StSbeaHer- 
mitian positive semi-definite (HPSD) matrix (note that 
it is Hermitian by construction) |19j . which results in a 
relatively cumbersome condition for T. The energy con- 
servation conditions for T can therefore be summarized 



(10) 



Lossless : T + = -2T^T, 
General : [-T - T* - 2Tf T] HPSD. 

III. THE ^-MATRIX 



A. Definition 

We here highlight an alternative formulation of the T- 
matrix method, which simplifies the energy conservation 
condition and naturally provides a connection with the 
radiative correction. Note that we will not here attempt 
to give a rigorous mathematical derivation, but rather 
focus on the new physical insights. Our proposed for- 
mulation is related to the reactance matrix or i-T-matrix, 
which can be formally defined as the Cayley transform 
of the ^-matrix |37| and has been previously discussed in 
the context of the general quantum theory of scattering 
[TH1 I3S"] . Explicitly, we have: 

K = i(I-S)(I + S)- 1 . (11) 

A simple consequence of this definition is that S unitary 
is equivalent to K Hermitian. In terms of the T-matrix, 
we have: 



K = -iT(I + T)" 1 = + T) -1 T. 
We note that T and K commute and we also have 



K + KT 



-iT = K + TK, 



(12) 



(13) 



which, to pursue the analogy with quantum scattering, 
may be viewed as the matrix version of Heitler's integral 
equations pfl [TH] [35] . T can be obtained from K using 

T = iK(I - iK)- 1 = i(I - iK^K, 



or from the following property: 



T 1 = iK 1 I. 



(14) 



(15) 



It is important to emphasize that the if-matrix and 
T-matrix formulations are fully equivalent from a formal 
point of view. However, in practice, since approximations 
are carried out in computing K and T (at the very least, 
truncation of these infinite matrices) , the equivalence is 
no longer strictly valid. We will in fact show that the K- 
matrix formulation is then the most appropriate one in 
approximate treatments where energy conservation needs 
to remain strictly enforced. This will lead us naturally 
to a generalization of the radiative correction procedure 
discussed earlier. 
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B. Energy conservation and the A'-matrix 

It is interesting to discuss the formal implications of the 
A'-matrix formulation for energy conservation. For non- 
absorbing scatterers, the unitarity of S, or Eq. [9]in terms 
of T, are equivalent to K being Hermitian: K = K^. 
For a ge nera l scatterer, the energy conservation condi- 
tion (Eq.[io|) can be shown to be equivalent to (iK^ — iK) 
being a Hermitian positive semi-definite matrix (for de- 
tails see App. [B]). More formally, this condition can be 
restated as K being a dissipative matrix [41] E2] • We can 
therefore rewrite the condition [TOl in terms of K as: 

Lossless : K = K^, 

General : K dissipative ([iK* - iK] HPSD) . 

These are much more natural conditions than the ones 
obtained for T (or for S). We note that (iKJ — iK) is 
simply, up to a factor of i/2, the skew-Hermitian part 
of K and the condition is therefore a generalization of 
Im(AT) > to the case where K is a matrix. We can 
therefore naturally identify the skew-Hermitian part of 
K as representing absorption while its Hermitian part is 
linked to scattering/dispersion. There is a clear analogy 
with simpler response functions such as the susceptibility 
X = e — 1 of a material or the polarizability a of a scat- 
terer, for which Im(a) corresponds to absorption (and is 
zero for lossless cases) and is subject to the condition: 
Im(a) > 0. In fact, the requirement that K be a dissi- 
pative matrix suggests that it is the matrix analogue of 
a scalar linear response function |8] like the polarizabil- 
ity a, and should therefore in addition satisfy causality 
and dispersion relations 43J akin to Kramers-Kronig re- 
lations. Such mathematical developments are however 
outside the scope of this work. We here only point out 
that in the case of non-absorbing scatterers, the condi- 
tions of optical reciprocity and energy conservation on K 
become closely linked (this is because symmetry and her- 
miticity become equivalent for real matrices, see App. [D] 
for further details). 



using the basis (M^, , JVl| 2) , N^ 2) ), where the latter 
two VSWFs use the (irregular) spherical Bessel functions 
of the second kind (which are superpositions of outgoing 
and ingoing spherical waves), in contrast to the usual 
spherical Hankcl functions of the first kind (which are 
outgoing spherical waves only). For this, we simply write 
Mf = Ml) 1 ' +«m£, , which separate the outgoing spher- 
ical wave VSWF into a sum of regular (M^) and irreg- 
ular (M' 2 )) contributions and obtain: 



E out (r) =E inc (r)+E sca (r) 
= + P,)M( 1 )(fc 1 r) + (b v + q^N^ihr) 

V 

+ i P M? ) {k 1 v) + iq u ^Xk 1 r). (18) 



The coefficients (b+q) m Eq. 17 can then be in- 
terpreted as the the sum of the incident field and the 
regular part of the scattered field; the latter can there- 
fore here be viewed as the regularized self-field, i.e. the 
non-diverging part of the field created by the scatterer at 
its own position (r = 0). The AT-matrix then represents 
(up to a factor) the linear connection between the expan- 
sion coefficients of the scattered field with those of the 
total field (incident + self-field). It can also be viewed 
mathematically as the linear connection between the ex- 
pansions coefficients of the irregular part of the outside 

("2) (2) 

field (i.e. those of M£, , N£, ; ) and those of its regular part 
(i.e. those of M^, n£^). This latter remark can be used 
to show that the A'-matrix can be computed as easily as 
the T-matrix in the most common implementation of the 
T-matrix approach, the EBCM [2TJI22]. Explicitly, we 
can show that (see App. [C] for details): 



K = PU 



(19) 



C. Relation to expansion coefficients 

We now discuss how the AT-matrix relates to the field 
expansion coefficients (in terms of VSWFs) . Recall that 
the T-matrix represents the linear connection (see Eq. [5]) 
between the field expansion coefficients of the scattered 
field E sca (in terms of outgoing spherical waves with 

and Nl> ) and those of the incident field 



VSWFs Mi 3) 

E 

N 

that 



lnc ^in terms of regular waves with VSWFs M 
(1) ). Multiplying Eq. 



(i) 



and 



13 



by the vector ( b ) , we deduce 



iK 



b + q 



(17) 



The physical meaning of this expression becomes appar- 
ent when we expand the total field outside the particle 



where we have introduced the matrix U such that Q = 
P + iU, which can be computed as easily as Q by sub- 
stituting spherical Hankel functions of the first kind, 
h n (x) = jn(x) + Wn{x), by irregular spherical Bessel 
functions y n (x) (note that j n {x) are the regular spheri- 
cal Bessel functions). In fact, iU can be viewed as the 
irregular part of Q, while P is its regular part [T5]. As a 
result, within the EBCM approach, K can be calculated 
as simply as T, if not more simply. 

One of the central themes of this work is to argue that 
the formulation of the EM scattering problem in terms 
of K is much more than a mere change of notation and 
presents in some cases many advantages in terms of both 
the practical implementations and the physical interpre- 
tations of the method. 
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IV. FORMAL DERIVATION OF THE 
RADIATIVE CORRECTION 

A. Radiative correction and self-field 

The observations in the last section provide a link with 
the radiative correction j3[ [5] , from the point of view of 
the self-field or self-reaction. Using again the example 
of a point polarizable dipole as illustration, the radia- 
tive correction can be interpreted as the effect of the 
self-field, i.e. the field Esf created by the scatterer onto 
itself, which acts in addition to the incident field E inc . 
The induced dipole is therefore written self-consistenly 
as p = ao(Ej nc + Egp), where ao again denotes the bare 
(uncorrected) polarizability. Since by linearity we have 
E SF = Gp (where G is the electric Green dyadic [J5] 
evaluated at the dipole position, taken isotropic for sim 
plicity), we obtain p = a 
polarizability a satisfies 



„RC 



Ei nc , where the corrected 



„RC 



) 1 =a 1 -G. 



(20) 



Classically, G diverges at the dipole position (which is 
why regularization is necessary [4]), but its imaginary 
part is finite and can be computed to recover Eq. [2j The 
if -matrix formulation provides a formal generalization 
of this approach. As mentioned earlier, the self-field is 
represented by the regular part of the scattered field, 
i.e. the part of its VSWF expansion including regular 

VSWFs, Ml 1) and only. Eq. 17 is the generaliza- 
tion of p = ao(Ei nc + Esf)- The formulation in terms of 
K therefore automatically includes this self-reaction. iK 
is analogous to Gao and represents the bare response, 
while T is analogous to Ga RC and corresponds to the 
self-reaction-corrected response, i.e. the radiative correc- 
tion. This analogy is further reinforced by comparing 
Eq. ~ 



15 and Eq. 20 



This crucial point can be further 
developed to provide a rigorous justification of the em- 
pirical radiative correction to the dipolar polarizability 
using the if-matrix, and by extending this concept to 
more general cases. 



B. Energy conservation in approximate treatments 

In practice, a number of approximations may be made 
when computing the T-matrix; at the very least, trun- 
cation is necessary. The energy conservation conditions 
(Eq. 10 1 may no longer be satisfied by the computed T- 



matrix, which is clearly undesirable (it could for example 
result in a predicted negative absorption cross-section). 
The optical reciprocity can be enforced a posteriori by 
appropriate symmetrization, but it is more difficult to 
enforce energy conservation. In contrast, the equivalent 
conditions on K (Eq. [Id] ) can more easily be enforced 
even when approximations are carried out (they are for 
example conserved upon truncation of the matrix). 



As an illustration, the T or K matrices may be approx- 
imated by expansions of their matrix elements [7], for 
example with respect to the size parameter (the lowest 
order being akin to the quasi-static or Rayleigh approx- 
imation), or with respect to the refractive index (more 
precisely, n — 1) for optically soft particles in the Born 
or Gans approximation [7j. In such instances, the energy 
condition on T may no longer be valid, notably because it 
mixes linear, e.g. T, and non-linear terms, e.g. T^T. On 
the other hand, the energy condition on K only contains 
linear terms in K and can be preserved. For example, the 
approximated K matrix will remain Hermitian for non- 
absorbing scatterers. We can therefore automatically en- 
force Eq. [lO] on T by deriving T from the approximated 
K using Eq. [14] or [Sj 

The reformulation of the scattering problem in terms of 
K therefore provides a natural method to enforce energy 
conservation in approximate treatments within the T- 
matrix framework, and this can for example be applied 
to the problem of radiative correction in EM scattering. 



C. General treatment of the radiative correction to 
the quasi-static approximations 

To illustrate further the procedure for deriving the ra- 
diative correction within the if-matrix approach, we fo- 
cus specifically on the important case of the radiative 
correction to the electrostatics and magnetostatics (or 
quasi-static) approximations. These approximations can 
be obtained from the general solution by taking the low- 
est non-zero order terms of the long-wavelength limit as 
k\ — > 0. We will denote T^ ) and K( ) the corresponding 
limit of the T- and Zf-matrices. These can in general be 
obtained from a direct solution of the electrostatics or 
magnetostatics problem, which is typically much easier 
than the full wave solution. 

As explained already, in general, the approximate T- 
matrix does not satisfy strictly the T-matrix energy 
conservation condition (Eq. 10); it only satisfies it ap- 
proximately to the accuracy to which it was calculated, 
i.e. in the long-wavelength limit as k\ — > 0. In contrast, it 
is relatively straightforward to ensure that the approxi- 
mate Zf-matrix K^ ^ satisfies exactly the X-matrix en- 
ergy conservation condition (Eq. 



16). If we therefore 
compute the approximate T-matrix from K^ - 1 , using for 



example Eq. 15 the resulting T-matrix will automati- 
cally satisfy the energy conservation condition and can be 
identified with the radiatively corrected T-matrix, T RC , 



(T RC ) _1 = -t (K<°) 



I. 



(21) 



These arguments provide a simple procedure to find 
the expression for the radiative correction for a given 
problem: 

• Solve the electrostatics and/or magnetostatics 
problem and find the corresponding "KS°\ which 



G 



should satisfy the energy conservation condition for 
K(Eq.[l6]). 



A. Physical interpretations of the vector spherical 
wave functions 



• Apply Eq. 21 to find the T-matrix with radiative 
correction. 

We note that for point scatterers, the first step is in 
fact implicit in the definition of its EM response, for ex- 
ample p = a E for an electric dipole. We also note that 
the matrix elements of T^ ' and K^°) are of order kf or 
higher. As a result, t he p roduct K^'T^ ) is of order at 
least kf, and from Eq. 13 we therefore have the following 



approximation for the matrix elements to lowest non-zero 
order: 



K (0) = -iT(°), 



(22) 



which can be used in many cases to find K' ) with stan- 
dard methods and results in: 



-.RCA 



T (0) 



This method can be used to generalize the concept of 
radiative correction to any type of scatterer, punctual or 
bodies of arbitrary shapes, including arbitrary multipole 
orders and interactions between multipolcs. Specific ex- 
amples that have been recently studied by other means 
include point multipoles (quadrupole, etc.) [H], bian- 
isotropic lossless point dipole scatterers [45 and point 
magnetic and electric dipoles with magneto-electric cou- 
pling [17]. The expressions obtained in all these studies 
are in fact special cases of Eq. 21 as we shall explicitly 
show in the next section. It is also interesting to note that 
Eqs. [15] and [21] can equally apply to study the radiative 
correction to higher order expansions of the polarizabil- 
ities in terms of k\ (as illustrated in the simple case of 
spheres in the next section) or to expansions in terms 
of other parameters (for example n — 1 for optically soft 
particles). 



V. APPLICATION TO SPECIFIC CASES 

We now study in more detail how the arguments pre- 
sented so far can be applied to specific cases of interest, 
some of which have been studied in the past using mostly 
heuristic arguments. As we shall see, all examples of 
radiative corrections studied so far in the literature are 

The only difficulty is to 



special cases of Eq. 15 or 21 



link the T-/A'-matrix formalism to more natural physi- 
cal representations in terms of for example polarizability, 
multipole moments, and multipole fields. We therefore 
first show explicitly that this link is relatively straight- 
forward; in essence it is simply a matter of definition and 
units. We then focus in the rest of this section on specific 
examples. 



The T- and Zf-matrix provide relations between coef- 
ficients of the expansions of the fields in vector spheri- 
cal wave functions (VSWFs). In practice, however, the 
excitation and response of the system are typically ex- 
pressed in a more natural form. For example, the exci- 
tation may be in the form of a constant external elec- 
tric field (in electrostatics) or the field of a plane wave. 
The response is often modeled in the form of an induced 
dipole (or multipole), the electromagnetic field of which 
implicitly represents the scattered field. For the applica- 
bility of the formalism, it is therefore necessary to link 
these physical excitations/reponses to their VSWF ex- 
pansions. The VSWFs of the scattered field, (M&m 

and 



N„m), are outgoing spherical waves and can be read- 
ily identified [3] with multipolar fields of order n (to- 
tal angular momentum) and angular momentum number 

(3) 

m. Nnm correspond to electric multipoles (also called 
(23) transverse-magnetic [3]), while M„m correspond to mag- 
netic multipoles (also called transverse-electric). The ex- 
pansion coefficients of the scattered field (p n m,qnm,) are 
therefore proportional to the magnetic and electric multi- 
pole moments of the scattered field (in a spherical tensor 
representation) . 

In a similar fashion, the expansion of the incident field 



in terms of a series of regular VSWFs, (M^N^l) can 
be viewed as the multipolar decomposition of the incident 
field, also in a spherical tensor representation. For an 
arbitrary incident plane wave, such an expansion can be 
computed analytically [19] and it is in fact one of the 
necessary steps in Mie theory [7J. 

In the context of the radiative correction to the electro- 
statics approximation, it is interesting to write explicitly 
these expressions. If we consider a general external elec- 
tric field (with sources at infinity), which is defined by 
an electric potential 0i nc (r) solution of Laplace equation, 
we may expand it as (the negative sign is chosen for con- 
venience) : 



:(r) = - 



iy 
1 n 



(24) 



where Y nm {9, (f>) are normalized scalar spherical harmon- 
ics (see App.[A|. The electrostatic response of a scatterer 
(point or body) to this external field can be written as a 
standard multipole expansion [3] of the potential created 
outside the scatterer as: 



0sca(r) 



1 



Y„ 



r n+l 



(25) 



The electric fields can be obtained from the standard re- 
lation E = — V0. In the general case, the induced mul- 
tipole moments (represented as a vector q) are linearly 
related to the excitation coefficients b by: 



ab. 



(26) 
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a is a generalized multipolar static polarizability tensor 
in the spherical tensor representation. We note that dif- 
ferent proportionality constants (potentially depending 
on n, m) could be introduced in the multipole expansions 
above and would affect the definition of a. In fact, the 
T-matrix formulation in the electrostatics limit (k\ — > 0) 
is an example of such an alternative definition. More 
explicitly, we can obtain the electrostatics limit for the 
normalized VSWFs (electric multipoles only) as: 



nm 



(2n + l)!!V n 



V(r"y nm ), 



N(3) 



i(2n-l)!! r~fi~ fY 7 



nm L.n+1 



V n + 1 \ r n+1 



(27) 



The electrostatics problem can therefore be recast within 
the T-matrix formulation as 



E- =Vi N (1) 

n,m 

E =Vo N (3) 

-'-'sea — / Hnm^nmi 



with 



f 22 b 



(28) 



(29) 



where T 22 is the electrostatics limit of T 22 (it is the bot- 
tom right block of T^). Note that T 11 , T 12 , and T 21 
involve magnetic multipoles and are zero in an electro- 
statics problem. 

Combining all this, we obtain a relation between the 
matrix elements of the multipolar static polarizability 
tensor and the electrostatics limit of the T-matrix: 



rpll 

nm 



ik 2n+l 



1 1 



(2n-l)!!(2rc+l)!! n 47re ei 



(30) 



The matrix T 22 is therefore simply, up to some propor- 
tionality factors, the multipolar static polarizability ten- 
sor. The radiative correction to the T-matrix (Eq. 21) 



therefore also applies to any definition of the multipo- 
lar polarizability tensors except for the proportionality 
constants (and correctly keeping track of these propor- 
tionality constants is the primary difficulty in writing it 
out explicitly). We will give specific examples in the fol- 
lowing. 

Finally, a similar result can be obtained for the mag- 
netostatics case in terms of the multipolar magnetic po- 
larizability tensor /3„ m , by substituting a nm = eoei/3 nm : 



c*o , typically obtained from a quasi-static (electrostatics) 
treatment. The response of such a dipole to an incident 
field Ej nc at its position is defined by the induced dipole 
moment: p = c*oEi nc . The T-matrix is here simply pro- 
portional to the polarizability tensor as given by Eq. [30] 
with n = 1: 



67reoe 



a, 



(32) 



which is valid both for the approximate T-matrix T^ ^ 
in terms of a.Q and the corrected T-matrix T RC in terms 
of a . The other blocks of the T-matrix are zero in 
this case. Moreover, T^ ) is of order k\ and we there- 
fore have K<°) = -iT^ (Eq. El. We note that for 



any physical polarizability tensor, c*o should be Hermi- 
tian if there is no absorption and a dissipative matrix if 
there is absorption (in a diagonal basis with eigenvalues 
CKj, this is equivalent to a% real, or Im(aj) > 0, respec- 
tively). K^ - 1 therefore satisfies the energy conservation 
conditions (Eq 



16 1. We can therefore rewrite Eq. 



21 



in terms of a to obtain the expression for the radiative 
correction as 



67re ei 



(33) 



which is the same as Eq. [2] previously obtained heuristi- 
cally for an isotropic polarizability tensor. The expres- 
sion above in fact extends it to the case of a general po- 
larizability tensor. Moreover, the argument remains valid 
for a body scatterer, when considering only the electric 
dipolar response. This for example justifies the empirical 
use of such a radiative correction for spheroidal particles 

nana. 

Note that in general cto depends on the frequency oj 
and the causality condition [S] for a RC cannot be eas- 
ily assessed by inspection of Eq. 33 only. We speculate 
that, based on the arguments in Ref. 43] the fact that K 
and its approximation are dissipative will automatically 
enforce causality for a , but further work (outside the 
scope of this paper) is necessary to investigate such as- 
pects. 

Finally, for a magnetic point dipole, whose response to 
an incident magnetic field Hj nc is an induced magnetic 
dipole moment m = /3 H inc , we would obtain following 
a similar reasoning: 



(/3 R T = 



(34) 



Jill __ 

Jim 



l l 



(2n- l)!!(2n + l)!! n 4tt 



(31) 



C. Electric point quadrupole and higher order 
multipoles 



B. Electric and magnetic point dipoles 

The simplest example of radiative correction is that 
of a point polarizable dipole with polarizability tensor 



It is straightforward to extend the arguments above to 
higher order multipoles. For electric multipoles, we can 
use Eq. |30| which relates the multipolar polarizability 
in the spherical tensor representation to the T-matrix. 
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Applying the general formula [2TJ we obtain the radiative 
correction for an electric multipole of order n as: 



i 4 U2n+1 

v ; v o; 47re ei (2n - l)!!(2n + 1)!! n 



h 

(35) 

which is the same expression as proposed in Ref. 03] 
for spheres. The expression above in fact extends it to 
the case of a general polarizability tensor (for example 
anisotropic). 



For magneto-electric coupling, the arguments of Sec. V A 
can be applied by noticing that to change from E to H 
one may make the substitution b — > (—i/Z)& for the 
incident field and q — > (—i/Z)p for the scattered field. 
We then get: 



ik\ 



-n=l 



rpl2 

A ri=l 



„SI _ 2 „-;J3, 



D7reoei Z o 

^(iZ)a s ^ H = ^ik\{icL HE ), (40) 



D. Point dipoles with magneto-optic coupling 

Sersic et al. have recently studied in detail the case of a 
magnetoelectric point-dipole [17 in the context of meta- 
materials and derived an expression for the radiative cor- 
rection (Eq. (18) in Ref. [T7]) using empirical arguments 
based on the optical theorem for this system [35]. This 
expression appears immediately similar to our general 
formula (Eq. 21) and the equivalence between the two 



can be demonstrated providing the prefactors and units 
are accounted for carefully, as we now show explicitly. 

Following Ref. [T7l the response of the scatterer to an 
incident electromagnetic field (E inc ,H inc ) is defined by 
the induced electric (p) and magnetic (m) dipole mo- 
ments obtained from the most general linear relation: 



Hi 



OLEE a EH 
OLRE OLHH 



(36) 



which is here written in rationalized units as in Ref. 1171 
a is a 6 x 6 polarizability tensor, compactly written in 
block-matrix notation. 

To apply our formalism, we rewrite this definition in 
SI units as follows: 



.SI 



m 



si 



nth 07] 



a SI = 47re ei 



si 



OLEE 



E SI 

nic 

ttSI 

ii- „ 



OLRE 



ZdLEH 



(37) 



(38) 



where Z = -\//io/( e o e i) = l/( e o c v^i) ^ s * rie impedance 
of the embedding medium (with relative dielectric con- 
stant £i). 

Using the arguments of Sec. V A[ we may express each 
block-matrix of the polarizability tensor in terms of a 
block-matrix (defined in Eq. [6]) of the T- matrix for n = 1 
(dipole terms). For example, using Eqs. 30 and 31 with 
n = 1, we obtain: 



rr\22 

± n=l 



ik\ 
67re ei 
ikf 

()7T 



EE 



-ik^EE, 



(39) 



These results can be written in more concise form as 
2 ., , ( a HH ioL HE 



T = 



lOLEH a-EE 



(41) 



Following our procedure presented in Sec. IV C we 
therefore obtain from Eq. [2l]the radiative correction for- 
mula for such a magneto-electric scatterer as: 



a 



RC 
HH 

.RC 
'■EH 



ICt 



RC \ 
HE 



„RC 
'EE 



a 



HH 



HE 



IOL 



EH 



a 



-ikfl. 



EE 



(42) 



To recast this expression in terms of the original defi- 
nition of a, it is necessary to change basis. One may 
introduce the unitary matrix W = ( j q ) and by left- 
multiplying by W and right-multiplying by W', we ob- 
tain: 



(a ) 1 - -ikfl, 



(43) 



which is the same expression as previously obtained em- 
pirically (Eq. (18) in Ref. IT?) . 



E. Beyond the electrostatics approximation 

As mentioned earlier, our procedure for radiative cor- 
rection can be applied to other types of approximations, 
for example beyond the electrostatic approximation. As 
an illustration, we will here study in detail the case of 
a spherical scatterer, for which exact results can be ob- 
tained from Mie theory [7j. In this case, the P-, Q-, U-, 
K-, and T-matrix are diagonal and independent of m, 
and the only non-zero terms are: 



pll 

nn 


= A n 


[sip n (x)ip' n (sx) 


~ ^'ni^nisx)] , 


^nn 


= A n 


[s£, n {x)lp' n (sx) 


~ £,n(x)i>n(sx)} , 


u 11 

nn 


= A n 


[sXn{xW n {sx) 


- Xni^nisx)} , 


pll 

nn 


= A n 


[lp n (x)ij' n (sx) - 


" S1p' n (x)tp n (sx)} , 


O 22 

^nn 


— A n 


[£,n{x)tp' n (sx) - 


s£,'n(x)Tpn{sx)] , 


u 22 

nn 


= A n 


[Xn{xW n {sx) - 




rpii 

nn 


pii 1 

r nnl ^nn > 




K ii 
nn 


pii IJjH 

nn/ ^ nm 





(44) 
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where 



n(n + 1) 



(45) 



x = k\a (with a the radius of the sphere), and s — 
'y/ei is the relative refractive index. The functions 
'4 ) n{x) 1 Xn{x) 1 and £ n (%) are the Riccati-Bessel functions 
[7] defined in terms of the spherical Bessel and Hankel 
functions as: 

1pn{x) = Xj n (x), Xn{x) = Xy n (x), 

Z n {x) = xh,W(x) = i) n {x) + ixn(x). (46) 

This is consistent with standard approaches to Mie the- 
ory since 



Pnm — fn^mj 
Qnm — A n 6 nm , 



(47) 



where T„ = = -P^JQW and A„ = T 22 = 

—Pn n /Q n 2 n are the magnetic and electric Mie suscepti- 
bilities: 



r„ 



sij) n {x)ijj' n {sx) - ijj' n (x)i> n (sx) 

s ^n{xW n {sx) - £' n (x)lp n (sx) ' 



A = ^n(x)j>' n (sx) - S^' n (x)lp n (sx) 
€n(x)lp'„( sx ) - s£' n (x)lj) n (sx) ' 

The extinction, scattering, and absorption coefficients 
(i.e. cross-sections normalized to the geometrical cross- 
section) can be obtained from these as: 



Qe 



r 2 Z-^i 



(2n+l)[Re(r„)+Re(A n )], 



^V(2n + 1) [|r„| 2 + |A„| 2 ], 



;sca n / . 

X z L — ' 



Qabs — Qcxt Qsca, ~ 

x [|r„| 2 Re(l + T- 1 ) + |A„| 2 Re(l + A" 1 )] . (49) 

Energy conservation, Q ext = Q sca + Qabs with <3 abs > 0, 
then requires that [40] : 



l + Re(A- 1 ) <0, 



(50) 



with the equality holding for non-absorbing spheres (for 
which s is real). Note that this condition is simply Eq. 
[T0| in the special case of spherical scatterers. In terms 
of the if-matrix, it takes a simpler form Im(i^"„) > 0, 
where K l n l n is up to a sign the same as r n and A„ upon 
substitution of $, n {x) by Xn{ x )- 

There have been many attempts to find suitable small- 
argument expansions of the Mie susceptibilities [5] EH 
|48l|49], notably in the context of plasmonics for the study 
of localized surface plasmon resonances (LSPR) in metal- 
lic nanospheres, where \s\ may be relatively large. The 
main dipolar LSPR is determined by Ai and its resonant 



character is evident in the wavelength-dependence of the 
far-field properties, which are then given by: 



G 



iext ~ -^Re(Ai), 



? sca ~ ^2 l Al l 



(51) 



The lowest order approximation to is 



(0)_2i s 2 -l 3 
Al ~ 3 s*+2 X ' [b2) 

which is simply equivalent to the electrostatics approxi- 
mation [7]. 

However, as illustrated in Fig.[T]for a silver nanosphere 
immersed in water, this approximation is only valid up 
to very small sizes of f=a 5 nm for metallic spheres, and 
in fact predicts a negative absorption as x increases. 
Moreover, since the electrostatics approximation is size- 
independent, it does not predict the redshift and broad- 
ening of the LSPR as the size increases. The radiative 
correction to this dipolar polarizability as given in Eq. [2] 
was in fact originally introduced empirically to remedy 
this problem [5]. It can be simply expressed as: 



Af 



(O)-RC 



= A? 



(0) 



- 1. 



(53) 



and is another example of an application of our general 
formula |21| However, as shown in Fig. [T] the improve- 
ment is marginal and not quantitative. It corrects the 
problem of negative absorption as expected, and predicts 
the strength of the resonance and its broadening, but not 
the size-induced redshift. 

Higher order expansions, up to third relative order, 
have been proposed, notably [4"9"] : 



Af = Ai' 



(0) 



1- 2-(.s 2 + l)+0(x 4 ) 



1 - 



10 s 2 +2 



i(s 2 + 10)- Ai 0) +0(x 



(54) 



In this expression, the numerator and denominator have 
been expanded to third order (relative to lowest order). 
As shown in Fig. [lj this significantly increases the range 
of validity of the approximation, up to a w 20 — 30 nm. 
For larger sizes however, although Eq. [54] predicts the 
correct redshift, it fails to predict the correct magnitude 
of the resonance. This can be attributed to the fact that 
A^ does not strictly satisfy the energy conservation con 
dition for the T-matrix (Eq 
here. 

The JC-matrix formalism here provides a simple 
method to address this issue and improve upon this ap- 
proximation. Instead of approximating directly the T- 
matrix (Ai = —Pii/Qii), we therefore use an approxi- 
mation of the if-matrix, which considering only the elec- 
tric dipole term is simply K\ = P 22 /C/ 22 . Expanding the 
numerator P 22 and denominator ?7 22 as before, we have: 



10), equivalent to Eq. 50 



K? = -iAi' 



(o) 



1- %(s 2 + l) + 0(x i ) 



i-!^(* 2 + io) + o(s*) 
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Lowest-order expansions Next-order expansions 

a - 5nm a = 10nm a = 20nm x ^ a = 30nm a = 50nm 

T ' 1 ' 1 1 1 I > 1 > 1 ' 1 ' 1 I ' 1 h 1 ' 1 ' 1 I — ' 1 ' 1 ' 1 ' 1 ' 1 I — ' 1 1 1 ' 1 ' r 




■ ■ ' ' 1 r | , i ■ I I ■ ■ t ■ ■ i I | , i , i , i "I I ■ ■ .1 

380 400 420 380 400 420 380 400 420 400 450 500 550 400 450 500 550 

Wavelength (nm) Wavelength (nm) Wavelength (nm) Wavelength (nm) Wavelength (nm) 



FIG. 1. (Color online) Predictions of the dipolar localized surface plasmon resonance for a silver nanosphere in water, as 
evidenced by the wavelength dependence of the far-field properties: extinction (Q C xt), scattering (Q S ca), and absorption (Q a bs = 
Qcxt — Qsca)- Only the dominant electric dipole response (corresponding to Ai) was included in these calculations. We compare 
the exact result (bold/blue lines) with approximate results for increasing sphere size. For the lowest sizes (radii of a = 5, 10, 
and 20 nm, we compare with the predictions of the electrostatics approximation (ESA) from Eq. 52 (red/dashed lines) and 
those of the radiative correction to the ESA fr om E q. |53| (green /solid lines). For larger sizes, we compare to the higher order 



51 



(pink/ dashed lines) and its proposed radiatively-corrected version A^ 



expansions approximations using A^ from Eq. 
from Eq. 56 (dark cyan/solid lines). Note that these higher order expansions are accurate for a < 20 nm (their predictions 
would lie on top of the exact results). In all cases, the vertical scale has been adjusted for best visualization of the quality of the 
approximation and the zero corresponds to the x-axis (except in the two cases where the hatched area indicates the negative 
region) . 



It is already apparent that the use of the fT-matrix pro- 
vides simpler expansions, as the denominator f/ff now 
has a well-defined parity as opposed to Qff (this is be- 
cause Xn(x) is odd or even, while t; n {x) is not). All odd- 
order terms in the denominator therefore disappear (they 
will in fact reappear as a result of the radiative correc- 
tion). Ki also satisfies the energy conservation condi- 
tion, Eq. [l6j at least for sufficiently small x. In fact, for 
s real (non-absorbing sphere), the condition is strictly 
satisfied for all x since is then real. We apply the 
central formula of this work, Eq. 21 to derive the radia- 
tive correction to the approximated and obtain: 



A-RC 



1 - 



10 s 2 +2 



(.s 2 + 10) 



(0) 



1 - 2M S 
1 10 V b 



1) 



- 1. 



(56) 



This expression differs from the earlier expression 



(Eq. 54 1 only by terms of relative order x 4 or larger. How- 



ever, because it was derived from the if-matrix formal- 
ism, it should be more physically valid, especially in res- 
onant systems where energy conservation is crucial. This 
is indeed the case, as it correctly predicts the LSP reso- 
nance behavior better than the previous approach based 
on a direct expansion of A^, in fact up to a « 50 nm as 
shown in Fig. [T] 



VI. CONCLUSION 

In conclusion, we have shown that the if-matrix pro- 
vides a simple formalism to study general radiative cor- 
rection problems in EM scattering. This was demon- 
strated in this work on two fronts. Firstly, from an ab- 
stract point of view, we studied and highlighted the for- 
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mal properties of the X -matrix, in particular with regard 
to energy conservation constraints. And secondly, from a 
practical point of view, we showed that previously pub- 
lished radiative correction formula are a straightforward 
consequence of the iVmatrix formalism when applying 
the method described in this work, namely: rather than 
obtaining approximations of the T-matrix, it is beneficial 
to derive it from an approximate IVmatrix using Eqs. |15| 
and 21 We expect that other systems will now be able 
to be studied following the same procedure. In addition, 
we believe that the iVmatrix formulation will play an 
important role in the general T-matrix approach to EM 
scattering. For example, as briefly mentioned in the text 
and in App. [Cj it provides an alternative route to nu- 
merical implementations of the T-matrix method, which 
may be more suited in some situations (for example for 
non-absorbing scatterers). 
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Appendix A: VSWF definitions 

The four types of VSWFs (j = 1,2,3,4) are defined 
using the same convention as Ref. [T^las follows: 

N(j) = J_v x M (i) 
ki 

(Note that : M<& = x N«), 
fei 



(Al) 



where 



7nr 



(2ti+1) (ri-m)! 



(A2) 



y Aim(n + 1) (n + m)\ 
is a normalization constant and 

(r, 9, <j>) = z^(kr) P™(cos(9)) e" n * (A3) 

are solutions of the scalar Helmholtz equation (with 

wave-vector k) in spherical coordinates, z^P are spheri- 
cal Bcssel functions, the choice of which defines the type 
of VSWFs (characterized by the superscript j): 

• Zn ^ — jn, i-e. spherical Bessel function of the first 
kind, for the regular VSWFs: j n {x) = x n /(2n + 

x— 

1)!! is regular at the origin. 

• Zn ^ — y n > i- e - spherical Bessel function of the sec- 



ond kind, for irregular VSWFs: y n (x) 
l)!!^; - " -1 is irregular at the origin. 



-(2n 



• z^ = hn — jn + iyn, he. spherical Hankel func- 
tion of the first kind, for outgoing spherical wave 
VSWFs: h ( n\x) ~ -(i) n+1 e ix /x. 

X — ^oo 

• z^ = hn — j n — iy n , i.e. spherical Hankel func- 
tion of the second kind, for ingoing spherical wave 
VSWFs: h ( n\x) ~ i n+1 e~ ix /x. 

x^-oo 

The associated Legendre functions P™(cos(6>)) are here 
defined with the Condon-Shortley phase, i.e. as: 

P™(x) = (-1H1 - x 2 ) m ' 2 —P n {x), (A4) 

where P n {x) are the Legendre polynomials. 

Note that in terms of normalized scalar spherical har- 
monics Y nm (9,4>) [3], we have: 



1 



V n ( n + 1) 



z^(k ir )Y nm (e^). (A5) 



Appendix B: Energy conservation condition for 
absorbing scatterers 

As mentioned in the main text, for absorbing parti- 
cles, the inequality cr cxt > cr sca requires that I — S^S 
be a Hermitian positive semi-definite matrix (HPSD) 
|19j . In terms of the T-matrix itself, this results in 
the somewhat cumbersome condition that the matrix 
- [T + Tt + 2TtT] is HPSD. 

A much simpler condition is obtained in terms of the 
iVmatrix by noticing that (assuming K and T are in- 
vertible): 



I - S f S HPSD 

& Kt [I - StS] K HPSD 

<S> -Kt [Tt(I + T) + (I + Tt)T]K HPSD 

O -KtTt(-iT) - (iTt)TK HPSD 

<S> Tt [iK+ - iK] T HPSD 

& [»Kt-*K] HPSD 



(Bl) 



The condition for the iVmatrix is therefore that 
[iKT — iK] is Hermitian positive semi-definite, or equiv- 
alcntly that if is a dissipative matrix [41] . 

Note that the same proof can easily be adapted to 
prove the special case that S unitary is equivalent to K 
hermitian. 



Appendix C: Computing the K- matrix with the 
Extended Boundary Condition Method (EBCM) 

One of the most common approaches to calculat- 
ing the T-matrix in practice is the Extended Bound- 
ary Condition Method (EBCM) or Null-Field Method 
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QH [HJ [H [Ml- Within this approach, T can be conve- 
niently obtained from T = — PQ _1 (Eq. rah, where the 
matrix elements of P and Q can be expressed analytically 
as surface integrals over the particle surface. Substitut- 



ing this into Eq. 13 and right multiplying by Q, we obtain 
K(Q — P) = iP . This leads us to introduce the matrix 
U such that Q = P 4- iXJ and we then have: 



K = PU 



(CI) 



In addition, the matrix elements of Q and P have iden- 
tical analytical expressions except for the substitution 
of M ( ^\k 1 r) for Q by M^far) for P. The matrix 
elements of U can therefore simply be obtained using 

(2) 

the same expressions but now with Mi (fcir) (this fol- 
lows from Q = P + tU and M^^r) = M^A^r) + 
iMffar)). Equivalents, U can be computed like Q, 
simply substituting spherical Hankel functions of the first 
kind, hn\x) — j n (x) + iy n (x), by irregular spherical 
Bessel functions y n (x). As a result, within the EBCM 
approach, K can be calculated as simply as T, if not 
more simply. 

The same conclusion could have been obtained by di- 
rect comparison of Eqs. [4][5] for T with the equivalent 
expressions [T7]|T8l 

In addition, it is interesting to note that the energy 
conservation condition for lossless scatterers is equivalent 
within the EBCM to 



U f P = P f U. 



(C2) 



Defining Y = U^P, this is also equivalent to Y Her- 
mitian. In the general case of absorbing scatterers, the 
condition that K be a dissipative matrix is also equiva- 
lent to Y dissipative. K can then be obtained from the 
following expression: 

PY -ipt 



K 



(C3) 



which automatically implies that K is Hermitian (dis- 
sipative) if Y is Hermitian (dissipative). These obser- 
vations allow one to check (and even enforce) the energy 
conservation condition on Y (for example as a function of 
truncation) before carrying out any matrix inversion. We 
believe that such an approach may also be further devel- 
oped to improve the numerical stability of the T-matrix 
approach, which is a common issue in EBCM implemen- 
tations 



reciprocity is equivalent to K 22 symmetric whereas en- 
ergy conservation is equivalent to K 22 Hermitian. For 
non-absorbing scatterers, the matrix elements of P 22 and 
U 22 are pure imaginary numbers and therefore K 22 is by 
construction a real matrix (from Eq. CI ). The conditions 
for optical reciprocity (K 22 symmetric) and energy con- 
servation (K 22 Hermitian) then become trivially equiva- 
lent. We note that this equivalence is not obvious when 
considering the T-matrix as opposed to the iC-matrix. 

This argument can in fact be generalized to the full 
if-matrix for a scatterer of arbitrary shape. The optical 
reciprocity condition then takes the form: 



n,m : n ,m 



= (-i) 



m+m' jy-ji 

n' ,—m' ,n. — rn ' 



(Dl) 



which is deduced from an identical relation for T (|19). 
Eq. 5.34). 

In the framework of the EBCM approach, K is com- 
puted from K = PU -1 , where the matrix elements of 
P and U are given by surface integrals involving cross- 
products of MW, NW, M( 2 ', M< 2 \ for example: 



Jn,rn,n' .m' ■ ^( l) 



dSn- 



M ( V , x M 

n'm' " 1 -n, — m 



r(D 



(D2) 

Moreover, in a non-absorbing medium (with wavevector 
k real), we have: 



M 



«_ m (fcr) = (-1)'" (M^(fcr) 



(D3) 



along with identical relations relating to N^, M^ 2 \ and 
N( 2 ) _ p or an integral like the one given above, we there- 
fore have: 



(J, 



n.m.n' ,m 



)* = -(-!) 



m+m 



n,— m,n' ,—m' 



(D4) 



By inspection of the integrals for the matrix elements 
of P and U (Ref. [TH] p. 145), we therefore deduce that: 



V n,m,r, 



= -(-i) 



m+m ptj 



n, — m.7i' ,—m' ' 



and 



n.m.n' ,m' 



m+m' jj'ij 



n, — m.n' , — m' 



(D5) 



(D6) 



By carrying out explicitly the block inversion of U and 
block-matrix multiplication of PU -1 , one may then show 
that 



Appendix D: Equivalence of optical reciprocity and 
energy conservation conditions on K 

We here restrict ourselves to non-absorbing scatterers. 
In order to highlight the central idea without being ham- 
pered by technicalities, it is enlightening to first consider 
the somewhat artificial case involving electric multipoles 
only, i.e. only the block T 22 of the T-matrix, and the case 
of a scatterer with symmetry of revolution (for which dif- 
ferent values of m are decoupled). In this case, optical 



(K ij , ,)* = (-1V 

V n,m,n' ,m' / V / 



n, — m,n' , — m' ' 



(D7) 



Using this expression (only valid for non-absorbing scat- 
terers), it is clear that the optical reciprocity condition 
(Eq. Dl) is equivalent to 



(K lJ , ,)* = K 3 ) , 

V n,7n,n,m / n,m,n,mi 



(D8) 



which is exactly the condition for energy conservation, 
K = Kt. 
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